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Abstract 

Nonadiabatic change of the control field or of the low-frequency coherence allows for an almost 
instantaneous change of the signal field propagating in a thick resonant absorber where electro- 
magnetically induced transparency is realized. This finding is applied for the storage and retrieval 
of the signal, for the creation of a signal copy and separation of this copy from the original pulse 
without its destruction. 
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A seminal idea of storage and retrieval of light pulses (SRLP) using a medium with 
electromagnetically induced transparency (EIT) was proposed in Ref. |l| by Fleischhauer 
and Lukin. Shortly after the proposal, SRLP was experimentally demonstrated in ultracold 
Na 2], in hot Rb vapor 3], and in a solid (Pr:YSO) [4j. Later, SRLP received much 
attention in the framework of quantum computing and quantum memory (see, for example, 
Ref. 5]). The core idea of such a storage is the dark-state polariton [lj. This is a particular 
superposition of the pulse amplitude and the atomic coherence created in a two-quantum 
process by signal and control fields. The signal pulse drives the transition from the populated 
ground state \g) to an excited state |e), and the cw control field drives the transition between 
state |e) and a metastable state |m), which are initially not populated (see inset in Fig. 
1). Reducing adiabatically the amplitude of the control field to zero, one can stop this 
polariton in a state that has zero amplitude for the signal pulse component and a nonzero 
amplitude for the component containing the coherence g — m. This coherence resembles a 
stand-still "spin wave" whose spatial shape coincides with the spatial shape of the signal 
pulse if it would have zero group velocity, V = 0. An adiabatic increase of the control 
field amplitude from zero back to its initial value retrieves the signal pulse from the "spin 
wave", both propagating with group velocity V ^ 0. Later, by numerical simulations and 
simple analytical calculations it was shown [6( that the adiabatic change of the control field 
amplitude is not crucial for SRLP. Even an abrupt change of the control field gives almost 
the same result: the stand-still spin wave is formed after the control field is switched off and 
then the signal pulse is retrieved after the control field is switched on again. 

Recently we showed [7] that actually the signal pulse entering an EIT medium is trans- 
formed at the very beginning to the control field and leaves the medium with group velocity 
c. Then an adiabaton is formed, which consists of a dip in the temporal profile of the 
control field and a bump for the signal field, propagating together with reduced group ve- 
locity V. For equal coupling constants and frequencies of the signal and control fields, the 
sum of their energies is constant in any cross section of the medium where the adiabaton is 
formed. Therefore, the bump and the dip in the intensity profiles of the propagating fields 
exactly compensate each other, and the energy of such an adiabaton is zero. This gives an 
important clue for the understanding of the slow light propagation. First, the spatial wave 
of the low frequency coherence (spin wave) is built up, propagating with the slow group 
velocity V. Then this wave commands the control field to produce the slow signal field, 
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transforming a part of the control field energy to and back from the signal field. Such a 
picture of the adiabaton suggests that, if we would instantaneously switch off the control 
field, the slow signal would not appear because no energy is available for the transformation 
described above. The spin wave produced at the early stage of the adiabaton formation is 
not destroyed but stops propagating because the source, i.e., the control field, which drives 
the spin-wave propagation, is switched off. Following this heuristic picture one can propose 
a new method to control a signal field. Firstly, one can instantaneously change the ampli- 
tude or phase of the control field. Since the control field is a source producing the slow 
signal field, its change would result in an instantaneous change of the amplitude or phase of 
the signal field. Secondly, one can almost instantaneously change the spin wave by a short 
rf pulse. Because the spin wave also controls the slow light production, its change would 
cause an instantaneous change of the signal field as well. This gives new opportunities to 
coherently control the properties of the radiation field. 

In this paper, we develop a novel scheme of such a control of the signal field and justify the 
heuristic arguments given above. In our scheme, first, we double the intensity of the signal 
pulse by the instantaneous doubling of the intensity of the control field. This operation also 
increases the group velocity of the signal pulse by a factor of two, i.e., it becomes 2V. Then 
by a short rf pulse we transfer half of the population of the metastable state \m) to another 
hyperfine level \M) of the ground state atom, also metastable. This imprints a snapshot 
of the spatial shape of the spin wave to state \M). Since state \M) is supposed to be not 
coupled to the other atomic states by the driving fields, this part of the spin wave becomes 
stand-still. But what is left in state \m) continues propagating with group velocity 2V. Its 
probability amplitude reduces by a factor of l/\/2 and, therefore, the intensity of the signal 
pulse drops by a factor of two back to its initial value. After some delay time, which is 
long enough to ensure that the signal pulse has already left the medium, we apply again an 
rf pulse to bring back the atomic population from state \M) to state |m), which has been 
already emptied by this time. The control field acts on the appearing stand-still spin-wave 
in state \m) such that the signal pulse is produced again and both, the spin- wave and the 
signal pulse, travel with the group velocity 2V. 

The excitation scheme is shown in the inset in Fig. 1. Initially, two ground state levels m 
and M are depopulated by optical pumping, and only state g is populated. The control field 
with the coupling amplitude (the Rabi frequency) O is a cw propagating along coordinate 
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z in a thick resonant sample. At time to = 0, the signal pulse with coupling amplitude \I/ 
enters the sample and propagates in the same direction as the field Q. It is assumed that 
\l/ -C f2 and the spectral width of the pulse A$ is smaller than the width of the transparency 
window, Ay = 2Q 2 / / y, where 7 is the decay rate of the coherence g — e, which is fast, 
7 > 2\/2Q. Since the signal pulse is weak, we can apply the linear response approximation 
for the solution of the Schrodinger equation for the atomic state 

|0) = C 9 \g) + C m \m) + C M \M) + C e \e) , (1) 

where Cm = 0. State \M) need only be considered when the rf pulse is applied. In this 
approach it is sufficient to consider only the evolution of the amplitudes C m and C e , which 
are described by the equations 

dX(z, t)/dt = t)Y(z, t), (2) 

dY(z, t)/dt = -^Y{z, t) - Q(z 7 t)X(z, t) - V(z, t), (3) 

where X = C m and Y = iC e . Here it is assumed that C g ~ 1 holds with a small deviation 
of the order of (^/Q) 2 1. If the condition of the adiabatic following of the dark state, 
A^ Ay, is satisfied, an approximate solution of Eqs. can ^ e easn y found [9( 

X(z,t) = -*(z,t)/n + ..., (4) 

Y(z,t) = -V t (z,t)/tt 2 + ..., (5) 

where ty t (z,t) = d^(z,t)/dt and the dots stand for terms that are at least A^/Ay times 
smaller. The wave equation for \1/ is 

L c V(z, t) = iaC* g C e w aY(z, t), (6) 

where a is a coupling constant and L c is the differential operator L c = d z + c~ 1 d t , where 
index c stands for the group velocity of the wave. Substitution of the solution (jSJ) into 
Eq. (jSJ) gives L c ^(z, t) = —(a/Q 2 )^/ t (z,t), which can be transformed to Lv 1 ^ r (z, t) = 0, 
where V\ = (c _1 + a/Q 2 )^ 1 is the new group velocity. The solution of this equation is 
*$?(z, t) = *Sf°(t — z/Vi), where ^°(t) is the amplitude of the signal field at the input of the 
sample. 

At time t\ > 0, when the pulse is in the sample, we abruptly change the amplitude of 
the control field: Q(z, t) = Q,[l+hQ(t — 1\ —z/c)], where 0(i) is the Heaviside step function 
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and we choose h = a/2 — 1 to double its intensity. This step-wise change of the control field 
amplitude propagates in the sample with velocity c. Before it arrives to the atoms with 
coordinate z, i.e., for t < t z = t\ + z/c, their amplitudes X(z,t) and Y(z,t) are described 
by Eqs. fll]),©. After t z , the solution of Eqs. ©,([3]) gives the following amplitudes 



Y{z,t) 



.'^A+K y (t-t z )^(t z -z/V 1 ), (8) 



where 



7 + e 7_T — 7_e 7+T 



K x (r) = 0(r), (9) 

7+ - 7- 



e" 



-7_r _ p-7+T 

Kir) = h - 0(r), (10) 

7+ - 7- 



7 ± = |±^-2fi 2 . (11) 

Here, the first and the second terms in the solution (El),® represent the main contribu- 
tions originating from the nonhomogeneous term, — ^(z, t), in Eq. and from the initial 
condition at t z . The omitted terms are at least A^/At times smaller. 

Before t z , the atomic state \<f>) was close to the dark state \d) = cos 8i(z, t) \g) — 
sm.9\(z,t)\m), uncoupled from the signal and control fields, where tan8i(z,t) = ty°/Q 
jlo| . The abrupt change of the control field amplitude from to y/2Q makes this state cou- 
pled since it acquires a particular component, which is the bright state \b) = sin 02(z, t) \g) + 
cos #2 (z, t) \m), coupled to the signal and control fields [lOj, where tan 9 2 (z, t) = *ff(z, t)/\/2Q. 
The probability amplitude of \b) is proportional to the signal field amplitude \1/ . Therefore, 
the amplitudes of the transient part of the solution, i.e., the second terms in Eqs. d7J),(jHD, 
are proportional to the signal field amplitude \1/ . Because of the fast decay of the excited 
state, the atom terminates its evolution in a new dark state with the mixing angle 02(z,t) 
[see Eq. ©]. 

With the solution (JSj) , the propagation equation (jBJ) can be transformed to 

L c ^(z, t) = -A{z, t)%{z, t) + aK y {t - t z )V°{t z - z/Vi), (12) 

where A(z,t) = fa + (fa ~ fa)®(t - t z ), fa :2 = V{£ - c^ 1 and V 2 = (cT 1 + a/21) 2 )- 1 is the 
new group velocity of the signal field after the jump of the amplitude Q(z,t). The solution 
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of Eq. ([12]) for t < t z is t) = - z/Vi). For t > t 2 this solution changes to 

#(M) = *°CO + <i>0M), (13) 

$(z,t) = a fXyit - t z , - Z -^-)^(t z , - z'/V x )dz', (14) 

JO V 2 

where T = «(V£/Vi)[t-t c - (z -z c )/V 2 ], k = {c-V x )/(c-V 2 ), and t z > = h + z'/c. z c = V 1 t c 
is a coordinate where at time t c = t\c/(c — Vi) the central part of the signal pulse changes 
its velocity from V\ to V 2 . Taking the integral (fl4|) by parts and retaining only the two main 
terms, we obtain 

$(z, t) = r]h^ (T) - V K x (t - t t )*°(t g - z/V x ), (15) 

where rj = (c — V%)jc. After a short time t — t z = r 7 ~ I/7, the function K x (t — t z ) decays to 
zero. If c ^> Vi, V2, then V2 = 2Vi, 77 « k « 1, z c « Viti and hence, for t > t z + r 7 we have 
t) = (1 + fc)* (T) where T = 2 [t - t x - - z c ) /V 2 ). This means that after the abrupt 
change of the amplitude of the control field the amplitude of the signal field also changes 
by the same factor (1 + h). In such a way the ratio ^(z,t)/Q(z,t) is conserved. Since 
C m (z,t) ~ —"$f(z,t)/Q,(z,t), the spin wave also conserves its amplitude and length. The 
latter coincides with the spatial length of the signal pulse in the sample before the change 
of the control field. This length is l p = V±t pi , where t pi = 1/A$ is the duration of the signal 
pulse at the input. Meanwhile, the spin wave and the signal field alter their group velocity 
from V\ to V%. Therefore, the duration of the signal pulse shortens to t P2 = t pl Vx/V 2 = t pi /2 
such that the spatial length, l p , of the pulse and the spin wave is conserved. 

At time t 2 = t\ + r 7 + l p j c, all spatial components of the spin wave and the signal pulse 
complete such a transformation. Following our scheme of the signal field processing, at 
time t 2 we apply a short, rectangular-shaped rf pulse, which drives resonantly the transition 
m — M (see inset in Fig. 1). The wavelength of the rf pulse is much greater than the spatial 
length, l p , of the signal pulse. Therefore, we disregard its spatial dependence. The evolution 
of the probability amplitudes C m and Cm of the atomic state \<p), Eq. (jTJ, is described by 
the equations 

dCJdt = iPC M + iVlC e , (16) 

dC M /dt = iPC m . (17) 

where P is the amplitude of the coupling m — M with the resonant rf pulse. We take 
P ^> Q and choose the duration of this pulse r r f such that it forms a so called it/2- 



pulse: Pr r f = ir/4 (see, for example, Ref. [llj for the definition). Before the rf pulse, 



we have Cm = and C m = — Since r r fQ <C 1, we can disregard the interaction 

with the control field during the rf pulse. Then, at the end of the pulse, £3 = £2 + T r /, 
we have C m = cos(Pr r/ )/fi = -#°/V2fi and C M = -z*° sin(Pr r/ )/fi = -i^°/V2Q. 
To simplify our consideration, we assume that the transition M — e is not allowed or far 
from resonance. Therefore the presence of the coherence m — M does not influence the 
signal and the control fields. Only the change of the probability amplitude C m introduces 
transients. They are described by Eqs. (J2J), (j3]), where *&(z, £) and Q(z, £) are replaced by 
y / 2 v I /0 and \/2fi (these amplitudes were present before the rf at t 2 ~ £3). At the end of the 
rf pulse, £3, the initial condition for an atom with coordinate z is X{z,t^) = — ^° (T3) / 
and Y(z,t 3 ) = -^° t {T)\ t=t ,jV2Q 2 , where T 3 = (VyVi)[£ 3 - *i - (« - O/V2]. After these 
modifications the solution of Eqs. (j2J), (EJ) is 

Y(z, £) « - ^(t - £ 3 )*°(T 3 ), (18) 

with the initial condition *t(z,i 3 ) = V2^(T)\ t=t3 . 

The solution of the wave equation (\6§ for £ < £ 3 is \&(,z, t) = ^/2^°(T). For £ > £ 3 , 
the r.h.s. of this equation changes to Eq. (fl8|) . Then its solution transforms to *&(z, t) = 
v / 2^r°(T') - $ 1 (2;,t) ) where $1(2, £) coincides with the function ®(z,t) in Eq. ([13]), if t z , 
and Vi are replaced by £3 and V2, respectively. Within the same approximation adopted for 
£), we obtain $i(z,t) = hr]^/°(T) — i]K x (t — t3)^°(T 3 ). Thus, after a short time r 7 , i.e., 
for £ > £ 3 + t 7 , and if 77 ~ 1, we have £) = \1/ (T). This means that after the switch off 
of the rf pulse the signal field changes its amplitude to its original value, which was at the 
input of the sample. 

We allow the signal field, having resumed its original amplitude, to leave the sample. The 
"spin wave", C m (z, £), accompanying the signal field and propagating with group velocity V2, 
vanishes at the sample end, z = l s . Meanwhile, the snapshot of the signal field at time £2 wa s 
imprinted by the rf pulse to the probability amplitude Cm(z,£ 3 ) = — i^/°(T 3 )/v / 2^ of state 
\M). This amplitude can be considered as another "spin wave" with zero group velocity 
(stand-still wave). Following the procedure described in the introduction, we apply a second 
rf pulse at time £4 when the signal field has already left the sample. The second rf pulse 
lasts six times longer than the first rf pulse, i.e., 6r r j, such that it forms a 37r-pulse: P6r r f = 
3%/2. The initial condition for this pulse is C m (z,t^) = and Cm{z, £4) = — i\l/ (T 3 )/v / 2fi. 
According to Eqs. f|T6|) . ffTTl) , at the end of this rf pulse, £5 = £4 + 6r r /, we have C m = 
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^°(T 3 )sin(P67; r/ )/v / 2fi = -^°(T 3 )/ v / 2fi and C M = -i^°{T 3 ) coa(P6r r f)/V2Q = 0. The 
extra 27r-rotation of the pseudospin 1/2, corresponding to the transition m — M, is necessary 
to obtain the proper sign for the final value of C m . Otherwise, if an rf pulse with 7r-area 
would be applied, the C m C* coherence would generate a signal field with a phase that is 
opposite to the initial one. In our case this coherence generates a field $ with the same 
phase as before. To show this we solve Eqs. with an arbitrary function $?(z, t) 

and tt(z,t) = y/2tt for the initial condition Y(z,t 5 ) = 0, X(z,t 5 ) = -^°(T 3 )/V2Q, and 
ty(z,ts) = 0. The solution for Y(z,t) is 

Substituting Y(z,t) to the wave equation ([6]), we obtain the solution 

*0M) = 1 1 K y {t-t 5 - —)^(T>)dz\ (20) 

where T 3 = (Vi/Vi)^ — t% — (z' — z^/V^. Approximately this integral is ^(z,t) = 
r)[5> (T 3 _ 5 ) e(t-t 5 )-*° (T 3 ) K x (t-t 5 )/h], where T 3 _ 5 = ^^[t-h+h-t,- (z-z c )/V 2 }. 
If 77 ~ 1, then for t > + r 7 we have \l/(z, t) = \1/ (T 3 _s) and the signal field is retrieved 
from the spin coherence. Now we have a copy of the signal field in the sample and the signal 
field outside the sample, both with the same amplitude and duration. Fig. 1 shows a 3-d 
plot of the signal field evolution controlled by the amplitude change of the control field and 
rf pulses. By numerical simulations we verified our approximate solution and obtained a fair 
agreement. It becomes almost perfect if our idealized solution is convoluted with a Gaussian 
function, described in 9j, which takes into account a pulse broadening due to the narrowing 
of the EIT window with distance. 

Summarizing, we found that the instantaneous change of the amplitude of the control 
field produces an almost instantaneous change of the amplitude of the signal field and it 
does not affect the amplitude of the low- frequency coherence (spin wave). This change also 
results in the variation of the group velocity and duration of the signal pulse. However, 
their product, which is the spatial length of the pulse and the spin wave, does not change 
if both are in the EIT sample when the change happens. The instantaneous change of the 
spin-wave amplitude by a short rf pulse produces an instantaneous change of the amplitude 
of the signal field without changing its group velocity and duration. By a train of rf pulses 
the signal field can be split in two parts, one of which can be temporarily stored in the 
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sample. These findings can be applied for information processing and storage, creation of a 
new type of entangled states, if a signal field contains only one photon. For example, a single 
photon can be split in two spatially and temporarily separated parts for the preparation of 
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time-bin qubits, which are of importance for quantum communication 
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FIG. 1: 3-d plot showing the spatial-temporal evolution of the signal field. Coordinates z, t and 
amplitude ^ of the signal field are in arbitrary units. The instances of the control, t±, t2, and £4, 
are indicated by arrows, as well as the end of the sample z = l s . The excitation scheme of the 
four-level atom is shown in the inset (see the text for details). 
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